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The general properties of the order parameter of a dipolar spinor Bose-Einstein condensate are 
discussed based on the symmetries of interactions. An initially spin-polarized dipolar condensate 
is shown to dynamically generate a non-singular vortex via spin-orbit interactions, a phenomenon 
reminiscent of the Einstein-de Haas effect in ferromagnets. 
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The realization of a Bose-Einstein condensate (BEC) 
of 52 Cr marks a major development in degenerate 

quantum gases in that the interparticle interaction via 
magnetic dipoles in this BEC is much larger than those in 
other spinor BECs of alkali atoms. The long-range nature 
and anisotropy of the dipolar interaction pose challeng- 
ing questions concerning the stability and superfluidity of 
the BEC @,i,H,@,0,lSll[Il3- The ground state of the 
52 Cr atom has a total electronic spin of three and a nu- 
clear spin of zero, and therefore the 52 Cr BEC has seven 
internal degrees of freedom. The interplay between dipo- 
lar and spinor interactions makes the order parameter of 
this system highly non-trivial [TlLIT3 .IT^. Moreover, the 
dipolc interaction couples the spin and orbital angular 
momenta so that an initial magnetization of the system 
causes the gas to rotate mechanically (Einstein-de Haas 
effect 0|) or J conversely, solid-body rotation of the sys- 
tem leads to its magnetization (Barnett effect [l^j)- 

This Letter investigates the Einstein-de Haas and Bar- 
nett effects in a spin-3 BEC system. We discuss the sym- 
metry of the order parameter of a dipolar spinor BEC 
and study the dynamic formation of spin textures using 
numerical simulations of the seven-component nonlocal 
mean-field theory, which takes into account short-range 
(van der Waals) interactions and magnetic dipole-dipole 
interactions subject to a trapping potential and an ex- 
ternal magnetic field. 

We first consider general properties of the order pa- 
rameter by discussing two fundamental symmetries of the 
dipolar interaction between the magnetic dipole moments 
fii = gfiBSi and fi 2 = g^BS 2 , where g is the electron g- 
factor, /ib is the Bohr magneton, and S\ and s 2 are the 
spin operators. The interaction between the magnetic 
dipoles located at r\ and r 2 is described by 



Vdd(ri 2 ) = c dd 



(si • s 2 ) - 3(si • e 12 )(s 2 • e 12 ) 



(1) 
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where r 12 = n - r 2 , e 12 = r 12 /r 12 , and c dd = 
/ i o(ffMB) 2 /47r with /zo being the magnetic permeability 
of the vacuum. The dipole interaction is invariant un- 
der simultaneous rotation in spin and coordinate spaces 
about an arbitrary axis, say the z-axis, so that the pro- 



jected total angular momentum operator S z + L z on that 
axis commutes with {i dd , where S = s\ + s 2 is the to- 
tal spin operator and L is the relative orbital angular 
momentum operator. Another symmetry of the dipo- 
lar interaction is the invariance under the transforma- 
tion P z e _l7rSz , where P z : (x,y,z) — ► (x, y, — z) and 
e~** S ' : (S x ,S y ,S z ) -> (-S x ,-S y ,S z ). Thus, the eigen- 
values of the following operators are conserved by the 
dipole interaction: 



S z + L z 



and 



-i-rrS z 



(2) 



A crucial observation is that these operators also com- 
mute with the short-range interactions. Thus if the con- 
fining potential is axisymmetric, the simultaneous eigen- 
functions of the two operators J5J) can serve to classify 
the two-body wave function. 

Constructing a many-body wave function by directly 
applying these symmetry considerations is quite compli- 
cated since the system has many degrees of freedom. 
However, substantial simplification can be achieved by 
considering the case of a ferromagnet in which dipole 
moments are localized at lattice sites and thus the de- 
grees of freedom of the system are reduced to the center- 
of-mass motion and the solid-body rotation around it. 
Consequently, spin relaxation of the system leads to a 
solid-body rotation of the ferromagnet — the Einstein- 
de Haas effect 14] . An analogous consideration can be 
applied to a BEC because almost all atoms occupy a 
single-particle state and therefore the degrees of freedom 
of the system can be represented by those of the order 
parameter. We may thus expect the Einstein-de Haas 
effect to emerge in a dipolar spinor BEC system. 

In general, the order parameter of a BEC can be de- 
fined as the eigenfunction corresponding to the macro- 
scopic eigenvalue of the reduced single-particle density 
operator. Let the order parameter be denoted as i/} a (r) 
with a norm assumed to be N, the number of condensate 
atoms. Here, a represents the magnetic sublevels of the 
atoms. It follows from the above symmetry considera- 
tions of the dipolar interaction that the order parameter 
of a dipolar spinor BEC in an axisymmetric system can 
be classified by the eigenvalues of the operators s z + l z 



2 



and P z exp(— iirs z ) corresponding to the operators J3J), 
where s z is the z component of the spin matrix and 
V = —i(d/d<f>). The simultaneous eigenstate of s z + V 
and P 2 exp(— ins z ) with eigenvalues J and p, respectively, 
is given by 



(3) 



where J is an integer and rj a j p is a complex eigenfunction 
of P 2 satisfying P z r/ aJp = p(-l) a r] aJp with p = ±1. A 
key point of Eq. © is that it includes the spin-dependent 
phase factor, since the dipolar interaction couples the 
spin with the orbital angular momentum. 

We now consider the dynamic formation of spin tex- 
tures in a dipolar BEC. The order parameter obeys the 
following set of the nonlocal Gross-Pitaevskii equations: 



ih 



d j> a {r) 
dt 



h 2 v 2 

' 2M 

2s 



gfi B B cxt a + C/trap(r) 



{ap\Vs\ot &)1>% (r) W (r)#< (r) 

Pa'P' 5=0 



E E B effW(^B^)#(r), 



(4) 



[i=x,y,z p 



where M is the atomic mass, B cxt is the external mag- 
netic field in the z direction, and t/trap^) is the trapping 
potential. We assume an optical trap so that all spin com- 
ponents experience the same trapping potential. The sec- 
ond line in Eq. (0J represents the short-range interaction, 
where the strength of the interaction is characterized by 
the s-wave scattering length a$ for a total spin S of a pair 
of atoms with spin s as g§ — Anitas /M . The operator 
Vs projects the wave function into the Hilbert space with 
a total spin S and is represented in terms of the Clebsch- 
Gordan coefficients {saaf3\SMg) as (a(3\Vs\oi' 1 ) = 
E S M s =~ S (sasP\SMs){SM s \ S a> S pi} 0. 

The last line in Eq. |@J represents the dipole-dipolc 
interaction, where s x,y,z are the spin matrices, and 



gm J \r 

E 

a'/3> 



„/|3 



rAr'K>0>^(r') (5) 



is the effective magnetic field at r produced by the 
surrounding magnetic dipoles, with e = (r — J"')/I r — 
r'\. Calculating the time derivative of S^(r) = 
Va( T ') S Q ( 3V'/3( r ')> we nn d that, apart from the spinor 
interactions, S(r) behaves like a classical spin and under- 
goes Larmor precession around the effective local mag- 
netic field B e s (r) -I- B cxt z. Hence, spin flip occurs in the 
region where B^(r) ^ 0. In a homogeneous infinite 
system, the effective field is completely canceled in a po- 
larized BEC and spin flip does not occur. Therefore, the 
Einstein-de Haas effect in an initially fully spin-polarized 
BEC is unique to non-uniform systems. 



We now examine the spin dynamics of the Einstein-de 
Haas effect in a spin-3 52 Cr BEC system. We consider 
a stable spin-polarized BEC in the lowest magnetic sub- 
level a — —3, produced in a strong magnetic field, as in 
the experiments of Refs. 0, Q We then suddenly de- 
crease the magnetic field to B ext . The initial state can 
be calculated by the imaginary-time propagation method 
in the subspace of ?/>_3(r). We have performed three- 
dimensional simulations of Eq. with seven spin com- 
ponents by using the Crank-Nicolson method. The scat- 
tering lengths of 52 Cr are reported to be = 112, 0,4 = 
58, and 02 = —7 in units of the Bohr radius [l7|. The 
value of ao is unknown and we estimate it using the 
van der Waals coefficient Cg given in Ref. ^t|> obtain- 
ing ao ~ (CqM/ttIc) 1 ' 4 — 91, where m c is the electron 
mass. Since we are interested in the spin dynamics of a 
fully spin-polarized state , th e short-range interaction is 
dominated by clq and 0,4 |lSj . 

We assume N = 10 5 atoms trapped in an axisymmetric 
potential U tTap (r) = (l/2)Mui 2 (r 2 + z 2 /\ 4 ) with ui = 
2ir x 820 Hz. The typical ratio of the dipolar interaction 
to the short-range interaction is s 2 Cdd/<?6 — 0-03 with 
s = 3. In a spherical trap, the number density at the 
trap center is n ~ 7 x 10 20 m~ 3 and the dipole-dipolc 
interaction energy s 2 Cddn becomes of the same order of 
magnitude as the Zeeman energy for B cxt ~ 0.1 mG. 
Hence, the spin-flip rate becomes significant for |B e xt| 

1 mG. In the following, we first consider B cxt = and A = 
1, and then discuss the effects of the external magnetic 
field and the trap geometry on the spin dynamics. 

Figure ^ shows the results for B oxt = and A = 1 . 
In Fig. OJa), we plot the population of each spin state 
N a /N = J dr\ip a \ 2 /N as a function of uxt. The figure 
shows that N-2/N first increases rapidly and then the 
components with a > —1 begin to increase. Figures^b), 
(c), and (d) show three-dimensional plots of 0_ 3 , -0-2, 
and ip—i, respectively, at uit = 2, where \4>a\ 2 is scaled 
by N/al a with Oho being the harmonic oscillator length 
y/H/2Mu). The order parameters show the symmetries of 
Eq. 0: t/j-2 has a phase factor and a node plane at 

2 = 0, and has a phase factor e -2 ^ and has reflection 
symmetry with respect to the z — plane. The other 
spin components have double-ring shapes similar to that 
of Fig. ^d) and their phase relationships satisfy Eq. © 
with J — — 3 and p = — 1. Beyond uit = 5, the spinor 
order parameter of the system develops a complicated 
structure. 

The effective magnetic field, Eq. at uit — and 
the spin vector S(r) at uit = 2 are plotted in Fig. The 
whirling patterns of the spin texture in Figs. |2fb) and 
(c) are due to Larmor precession around the local mag- 
netic field, shown in Fig.^a). Since the local magnetic 
field points outward for z > and inward for z < 0, 
the directions of the whirlpools in Figs.EIb) and (c) are 
opposite. Topological spin textures in spinor BECs have 
been observed in a spin- 1/2 system [l9| and in a spin-1 



3 



(a) 



(b) 



0.6 



0.4 



~\ 1 


1 


a=-3 

-2 


4.1 






-1 

o " 








1 

2 

3 










1° 




i-'-":" " -' : 




N 




FIG. 1: (Color) (a) Relative population N a /N = 
J dr \ip a \ 2 /N of each magnetic sublevel a versus ut for B ey ± = 
and A = 1. (b)-(d) Isopycnic surfaces of (b) tps, (c) ip-i, 
and (d) V-i at ut = 2, where \ip a \ 2 ai /N = 0.0001 for (b) 
and (c) and 5 x 10" 5 for (d). The color on the surfaces rep- 
resents the phase of the order parameter (refer to scale at 
right). 



system [2(j. We note that the spin textures in Figs.|2Jb) 
and (c) are generated spontaneously due to the intrinsic 
dipole interaction, in contrast to those in Refs. (lil l20|. 
which are generated by external forces. 

When the external magnetic field -E> xt is applied in the 
positive z direction and is much stronger than the dipole 
field, the spin angular momentum should be conserved 
because of energy conservation and spin flipping is sup- 
pressed. When B cxt < 0, spin flip can occur by convert- 
ing Zeeman energy to kinetic energy. These behaviors 
are demonstrated in Fig.^a), which shows the time evo- 
lution of N- 2 /N for B cxt = ±1 mG. Figures^) and (c) 
show cross sections of |-0_2| 2 af lo /^ for B cxt — ±1 mG 
at ut — 8 and A = 1. When £> GXt > 0, ip-2 oscillates in 
time between the structures of Fig. ^c) and Fig. Otb). 
These structures derive from the symmetry of the dipole 
interaction, which can be expressed in terms of rank-2 
spherical harmonics Yi m . The dipole field produced by 
an approximately spherical distribution of i/j_3 ~ Yqo(v) 
induces a Y2-i(r) term in ip-2 (Fig. Hfc)), which in turn 
affects itself and induces a linear combination of l^_i(r) 
and yj_i(r), resulting in Fig.|3Jb). Therefore, the struc- 
ture in Fig. [3ft)) manifests as a secondary effect of the 
dipole-dipole interaction. In the case of B cxt = — 1 mG, a 
similar structure as in Fig.^b) appears for ut < 4. How- 
ever, as time advances, the domain structure develops as 
shown in Fig. Etc)- The domain size becomes smaller as 
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FIG. 2: (Color) (a) Dipole field at t = for A = 1. The solid 
line shows the isopycnic curve at |-i/>_3| 2 a| /iV = 0.0001. The 
color of the arrows denotes the magnitude of the field. (b),(c) 
Spin configurations on the (b) z — 2 /im plane and (c) z = 
— 2 /im plane at ut = 2 for _B cxt = and A = 0. The length 
of the arrows represents the magnitude of the spin vector 
projected on the xy-plane and the color shows \S\a^ /N. Note 
that the spins tilt in a direction perpendicular to B c g(r). 
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FIG. 3: (Color) (a) Population of a = —2 sublevel versus uit 
for _B ex t = ±1 mG and A = 1 and 0.5. (b),(c) Cross sections 
of \ip- 2 \ 2 al /N with A = 1 at ut = 8 for (b) B cxt = 1 mG and 
(c) _B cxt = — 1 mG. The cross sections are at z = 0.7 /im (left) 
and y = /im (right). The size of each panel is 8.2 x 8.2 /xm. 



| B xt | increases. 

Finally, we discuss the geometry dependence of the dy- 
namics. Figure 0] shows the dipole field B c g(r) at t = 
for A = 0.5. Compared with Fig. Ufa), the z compo- 
nent of the effective magnetic field is inverted around the 
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FIG. 4: (Color) Dipole field at t = for A = 0.5. The solid 
line shows the isopycnic curve at \tp-3\ 2 a^ /N = 0.0005. Note 
that the sign of B^ a (r) in the condensate is opposite to that 
in Fig. 2(a). 

center of the condensate. The transfer of atoms from 
the spin component a — —3 to a > —2 is due to Lar- 
mor precession caused by B^, and it occurs most effi- 
ciently where the local magnetic field in the z direction 
B*g + Bext vanishes. The sign of such an optimized B cxt 
for spherical traps is opposite to that for pancake-shaped 
traps, as can be inferred from Fig. Ufa) and Fig. 0] This 
fact is reflected in Fig. Ofa) as the difference in the field 
dependence of the initial peaks. Since the z component 
of the dipole field is positive for most of the condensate 
when A = 1 (Fig. the initial peak in Fig. EJa) is 

larger for B cxt = — 1 mG than for B cxt = 1 mG. The 
relation between the initial peak and B cxt is opposite 
for A = 0.5, since the z component of the dipole field is 
mostly negative in the condensate (Fig. . In the case of 
a cigar-shaped BEC, the qualitative behavior is the same 
as for the spherical trap. 

In conclusion, we have shown that the Einstein-de 
Haas effect occurs in dipolar spinor Bose-Einstein con- 
densates and that a non-singular vortex appears from 
an initially spin-polarized condensate. In a low magnetic 
field (~ 1 mG) the fraction of spin-flipped atoms (~ 5 %) 
is large enough to be observed in Stern-Gerlach exper- 
iments. The spin-relaxation processes produce various 
vortex structures, depending on the external magnetic 
field. 
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Note added. — Very recently, a preprint |2l| has ap- 
peared which also discusses the Einstein-de Haas effect 
in a dipolar BEC. 
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